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1. Introduction 



In nuclear physics, the symmetry energy is first encountered within the 
elementary Bethe-Weizsacker formula for nuclear energy: 

Z 2 (N — Z) 2 

g; E = -a v A + a s A 2 ^ + acj^+a A - ^~^+ A 5 M 

it is the change in nuclear energy associated with changing neutron-proton 
, asymmetry (N — Z)/A. In nuclear matter, the energy per nucleon, depen- 

\ dent on neutron p n and proton p p densities, may be represented as a sum 

of the energy Eq for symmetric nuclear matter and the correction E\ asso- 
ciated with the asymmetry, 

E(Pn,P P ) = E a (p) + E 1 (p n , p p ) , (2) 

where p = p n + p p . The charge symmetry of nuclear interactions, which is 
the symmetry under the interchange of neutrons and protons, requires that 
the correction be quadratic in the asymmetry, for small asymmetries: 

El =E-E^S(p){^^>j 2 . (3) 

Microscopic calculations, such as, 1 indicate that the quadratic approxima- 
tion yields a very good representation for the energy in nuclear matter, up 
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to the limit of neutron matter with asymmetry of 1, over a wide range of 
net densities. As a consequence, the energy in nuclear matter over a broad 
range of parameters can be described exclusively in terms of the energy of 
symmetric matter Eo(p) and the symmetry-energy coefficient S(p). 

Much effort has been dedicated in the past to the understanding of 
Eo(p) and less so to S(p) whose features remain more obscure. 2 The en- 
ergy Eq(p) minimizes at the normal density po, reaching there the value 
of Eq w —16.0 MeV. The uncertainties in S hamper predictions for neu- 
tron stars whose structure depends on pressure in neutron matter. 3 In 
pure neutron matter, the energy is E = Eq + S and the pressure is 
P = p 2 dE I dp ~ p 2 dS/dp close to po, as Eq minimizes at po. In the calcu- 
lations of neutron-star structure, 3 a correlation is found, RP^ 1 ^ w const, 
between the radius R of a neutron star of a given mass and the pressure P 
in neutron matter at a given density p ~ pq. 



2. Binding Formula 

When examining the standard energy formula Q , we notice that the sym- 
metry energy has a volume character: it changes as A when the neutron and 
proton numbers are scaled by one factor. The formula lacks a surface sym- 
metry term that would change as A 2 ^ 3 with the change in nucleon number. 
A question to ask is whether there should be such a term. Let us look at 
the surface energy. We can write it as E$ = as A 2 / 3 = 4tt r 2 A 2 ! 3 — 

■ gS 2 S. The ratio ^M- = a = A as a is surface tension, the work that needs to 

47T 7-q 6 47T r£ ' 

be done per unit area when changing the surface area of the nucleus, such 
as in deforming the nucleus. The work needs to be done to compensate 
lost binding, as nucleons close to the surface are less bound, due to fewer 
neighbors, than in the interior. The formula states that in the interior 
the nucleons are less bound in a more asymmetric nucleus. In that case, the 
energetic price for increasing the surface should drop, i.e. a = ^£ should 
decrease with asymmetry, in the more general definition of tension. 

As intensive, a should be expressed in terms an intensive quantity as- 
sociated with the asymmetry, which is the asymmetry chemical potential: 

» A= d(N-Z) = 2 { ^-^- (4) 

To lowest order, under charge symmetry, the tension needs to be quadratic 
in p A , 

(T = a -jp 2 A . (5) 
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If the tension depends on asymmetry, so must surface energy. On examining 
the function $ = ^a{N — Z) — E, with the derivative d^/d^A = N — Z, the 
last dependence is seen to produce the following apparent paradox: some 
of the nuclear asymmetry N — Z must be associated with the surface and 
not the interior. To answer the question on how particles can be attributed 
to the surface, one needs to adapt a systematic approach the separation 
of quantities into volume and surface contributions. Gibbs 4 proposed to 
consider two copies of the system, actual and an idealized reference copy 
where the interior densities of different quantities extend up to the surface 
position, cf. Fig. ^ The idealized system represents one with only volume 




r 



Idealized 



R r 

Fig. 1. Gibbs' 4 construction for defining volume and surface contributions. 



contributions to different quantities, while the difference between the sys- 
tems can be associated with the surface, F$ = F — Fy. The separation, 
however, depends on the position of the surface which must be set utilizing 
some auxiliary condition. For nuclei, it is natural to demand a vanishing 
surface nucleon surface number As — 0, i.e. set the surface position at the 
sharp-edged sphere radius R. Since nuclei, though, are binary systems, the 
surface positions that might be separately attributed to neutrons and pro- 
tons will be, generally, displaced relative to each other, cf. Fig. O In conse- 
quence, in spite of a vanishing surface nucleon number, As = Ns + Zs = 0, 
we may have a finite surface asymmetry Ns — Zs ^ 0. 
Having resolved the apparent paradox, from JSJ we find 

77 c i 2 c Z70 i 1 ( N S - Z s) 2 P o I S ( N S - Z S ) 2 

E s =a Q S + ~iii\S = E s + — - =E s +a A . (6) 

Here, the index refers to symmetric matter and we have introduced a 
surface symmetry coefficient a^, with dimension of energy. For the volume 
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Fig. 2. In a two-component system, the surfaces for the two components will be, 
generally, displaced relative to each other. 



energy, within the standard formula, we have 



Ey = Ey 



.4 



(7) 



where a\ is the volume symmetry energy coefficient. The Coulomb energy 
is temporarily ignored. The net energy and asymmetry are, respectively, 
E = E s + E v and N - Z = N s - Z s + N v - Z v . In the ground state, 
the asymmetry should partition itself into the surface and volume contri- 
butions, in such a way as to minimize the energy. The result of the energy 
minimization can be, actually, written right away once one notices that the 
surface and volume energies, quadratic in asymmetry, are analogous to the 
energies of capacitors quadratic in the charge. The energy of the coupled 
capacitors is quadratic in the net charge, with the square divided by the 
net capacitance, yielding: 

- 2 (JV - Zf 



E = E U 



^ = E° 
2C 



A 2/3 



(8) 



Adding now the Coulomb term, we arrive at the modified energy formula 



E 



-a v A + a s A 2 ' 3 + a c — — + 



A i/3 i + A-VSoX/a 



(N -Zf 
A 



0) 



Compared to the standard formula lJT)l. the symmetry coefficient becomes 
now mass dependent, cla{A) — a\/{l + a\/(a^A 1 ^ 3 )). The standard for- 



mula is recovered when a\ 



/a* 



A ~ 0, i.e. the surface does not accept asym- 
metry excess, or in the limit of A — > oo. Within the modified formula, the 
symmetry coefficient weakens at low A. 5,6 Whether or not the coefficient 
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may be replaced by the constant a\ depends on the ratio a\/a^. That 
ratio can be determined from the ratio of surface to volume asymmetry 
partitioning for the energy minimum in proportion to the capacitances: 

N s -Zs _C S _A^/a s A ._,„a v 



= A- 1 '^. (10) 



N v -Z v C v A/a v A a s A 
3. Asymmetry Skins 

Establishing the relatively small differences in the distribution of neutrons 
and protons in nuclei has been, generally, difficult experimentally. Probes 
with different sensitivities to protons and neutrons had been utilized, such 
as electrons and protons, negative and positive pions, or protons and neu- 
trons, with different associated systematic errors. The results have not been 
expressed in terms of the surface excess, but rather in terms of the differ- 
ence in the r.m.s. radii between neutrons and protons. The conversion from 
the excess (|10f) to the difference of radii is relatively straightforward, 6 if 
the surface diffuseness is similar for neutrons and protons. Another issue 
to consider theoretically is that, for heavy nuclei, Coulomb forces compete 
with symmetry-energy effects, pushing the proton radius out against neu- 
trons and polarizing the nuclear interior. That competition is easily taken 
into account by minimizing the sum of three energies with respect to the 
asymmetry: 



E = E V + E S + E C , (11) 



where 



Zc = ~[^Z* + Z v Z s+ 1 -zi). (12) 

From the energy minimization, an analytic formula for the difference of the 
radii follows, 

(r% /2 - (r 2 )^ 2 A N-Z 



(r 2 )V2 QjSfZ 1 + A 1 / 3 a A /a\ 

ac A 5/3 m +A l/3 a S A/a V 

168< N l + i4Va a ^/oX 



(13) 



The first term on the r.h.s. represents the effects of symmetry energy only, 
from Eq. (|l(Jfl . while the second term represents the Coulomb correction. 
It should be mentioned that the impact of the Coulomb-symmetry energy 
competition is much weaker onto the net energy than onto the skin size. 
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Before trying to draw conclusions from data with Eq. it may be 

worthwhile to test the macroscopic theory against the microscopic. In their 
nonrelativistic Hartree-Fock and relativistic Hartree calculations, Typel and 
Brown 7 observed correlations between the sizes of asymmetry skins in dif- 
ferent nuclei, when utilizing different effective interactions. Those correla- 
tions are shown in Fig. |3| together with the predictions of Eq. Ill -ill when 
changing the ratio a~y/a\. The accuracy of the macroscopic theory in re- 
producing correlations from the microscopic theory appears to be at the 
level of 0.01 fm! 




/ 2\l/2 / 2\l/2 It \ t 132 i 138 n 

(r ) n — \r ) p (fm) for Sn and Ba 

Fig. 3. Correlation between the asymmetry skins for 208 Pb and 132 Sn and 138 Ba, in 
the nonrelativistic and relativistic mean-field calculations 7 (symbols) and predicted by 
macroscopic Eq. 1131 (lines). 



We next turn to the implications of skin data. Figure 0] shows a com- 
parison of the data by Suzuki et al. from Na isotopes to the predictions of 
Eq. (jT3|) for different values of a^/a\. The comparison suggests a value of 
a l/ a A ~ 3- Figure El shows, with sloped parallel lines, the constraints on 
a^/a\ from fitting a variety of skin data (for references to the experiments 
see 6 ). Without Coulomb effects the constraint lines would have been hori- 
zontal; the weak sensitivity of the fits to a\ results from the second term 
on the r.h.s. of (|13|l . The favored values of the ratio, a^/a\ ~ 2.8, imply 
that A~ 1 I Z a\/ a\ is never small. Neither can the A-dependent symmetry 
coefficient be replaced by a^, nor even expanded linearly in A -1 / 3 . 
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Fig. 5. Constraints on the symmetry energy parameters in the plane of a s A /a}^ vs a^. 
The sloped horizontal lines represent constraints on aj| / from fitting skin data at an 
assumed value of aX A . The elliptical contours represent constraints obtained from the fit 
linear in A -1 / 3 to the values of (A) from IAS. 
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4. Isobaric Analogue States 

To find absolute values of the symmetry coefficients, one might try to fit the 
binding formula to measured energies. However, this is treacherous as con- 
clusions on details in different isospin-dependent terms, including Coulomb, 
Wigner and pairing get interrelated when drawn from a global fit to the 
energies. In addition, the correlation between mass number and asymmetry 
along the line of stability correlates the conclusions on details in the isospin 
dependent and isospin independent terms. The conclusions on the symme- 
try coefficients change depending on what is done to the other terms in the 
formula. 6 

Optimal for determining the symmetry parameters would be a study of 
the symmetry term in the binding formula in isolation from the formula 
remainder, which might seem impossible. However, one can take advantage 
of the extension of charge symmetry of nuclear interactions to charge in- 
variance. Under charge invariance the symmetry term should be a scalar in 
isospin space 9 and can be, thus, generalized with 

E A = a A (A) ( N -Z? _ 4 a A (A) S 

T 2 T(T + 1) 

-^4a A (A) — = 4 a A (A) { ^ ) , (14) 

where we also happen to absorb most of the Wigner term into the symmetry 
term. Under the generalization, the binding formula may be applied to the 
lowest state of a given isospin T in a nucleus. When excited, such a state 
is an isobaric analogue state (IAS) of the ground state of a neighboring 
nucleus. In the formula generalization, the pairing contribution depends on 
the evenness of T. 

For an isospin of the same evenness as the ground, the change in the 
formula in the excitation occurs only in the symmetry term: 

E 2 (T 2 ) - E 1 (T 1 ) = ^{T 2 (T 2 + 1) - + 1)} , (15) 

and the excitation energy can be used to the determine the symmetry en- 
ergy nucleus by nucleus from 

AAE 

In the context of the previous considerations, the question to ask is whether 
the deduced ^4-dependent symmetry coefficient weakens for light nuclei and 
whether the inverse of the coefficient is linear in A -1 / 3 : 

a A \A)l{a v A )- 1 + {a s A )- 1 A- 1 /\ (17) 
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Fig. 6. Inverse of the A-dependent symmetry coefficient as a function of A -1 / 3 . Circles 
represent values extracted with 1151 from extremal IAS excitation energies in. 10 Circle 
size is proportional to the factor AT 2 /A in the coefficient determination. The line and 
squares show results of the fits to the experimental results following either Eq. 1171 or 
Thomas-Fermi theory. 13 

Inverse values of symmetry coefficients, extracted according to Eq. i|15fl 
from IAS data, 10 are shown in Fig. |S| It is seen that the inverse coeffi- 
cient changes with A -1 / 3 in a roughly linear fashion, although significant 
shell effects are present. The line across the figure represents best fit with 
Eq. (|17|l . The 1- and 2-a constraints on the symmetry coefficients from the 
fit are further indicated with elliptical contours in Fig. [5j Combining the 
constraints from the fits, we conclude that 30.0 MeV < a\ < 32.5 MeV and 



5. Consequences and Conclusion 

The emergence of the surface capacitance for asymmetry may be tied to 
the weakening of the symmetry energy with density, see e.g. 6,11 Due to the 
weakening, it becomes advantageous for the nucleus to push its asymme- 
try to the surface to lower energy. The ratio of the symmetry coefficients, 
specifically, can be tied to the shape of the symmetry energy dependence 
on density as, in the local-density approximation to the symmetry energy, 
the ratio is found to be 



2.6 < a\/a s A < 3.0. 




(18) 
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Here, the integration is across the nuclear surface and p(r) is the density as 
a function of position. For density- independent symmetry energy, S(p) = 
S(po) = a^, the surface does not accept the asymmetry, a^/a^ = 0! 

Using the correlations between the coefficient ratio, skins and drop of 
the symmetry energy with p, within the relativistic and nonrelativistic cal- 
culations by Fuhrnstahl, 12 one can arrive at limits at on the drop, either 
expressed in terms of the value of symmetry energy at half of the nor- 
mal density or in terms of the power of density in parameterization of the 
symmetry energy. 13 Specifically, one finds 0.58 < S(po/2)/a A < 0.69 and 
0-54 < 7 < 0.77 in S(p) ~ a\(p/po) J . These further imply limits on pres- 
sure in neutron matter at normal density and, with results of Ref., 3 produce 
limits on neutron-star radius of 11.5 km < R < 13.5 km for 1.4 Mq mass. 

To conclude, the requirement of macroscopic consistency brings in the 
surface symmetry energy into the nuclear binding formula. The volume and 
surface symmetry energies combine as energies of coupled capacitors. The 
extension of the binding formula implies emergence of the asymmetry skins 
for nuclei and weakening of the symmetry term in light nuclei. The sys- 
tematic of the asymmetry skins restricts ratio of the symmetry coefficients. 
The charge invariance allows to study variation of the symmetry coeffi- 
cient nucleus by nucleus. Combination of the fits to skins and IAS yields 
30.0 MeV < a\ < 32.5 MeV and 2.6 < a v A /a s A < 3.0. The surface sym- 
metry energy is associated with weakening of the symmetry energy with 
density. The a\/a^ ratio implies limits on drop characteristics, such as 
0.58 < S(po/2)/a A < 0.69. Implications for neutron stars follow. Current 
direction is to incorporate shell corrections into the IAS analysis. 
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